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Abstract: Using geometric invariant theory and momentum map geometry we show 
how to effectively find all stochastic local operations and classical communication 
(SLOCC) classes of pure states for both distinguishable and indistinguishable parti- 
cles. We prove that they are parametrized by critical sets of the total variance of the 
state which can be regarded as an entanglement measure. Remarkably, the Morse in- 
dex at the critical point have a nice interpretation of number of non-local directions in 
which such defined entanglement increases. We also introduce the SLOCC-invariant 
'measure' of entanglement as square root of total variance of state at the critical point 
and show how it is connected to the distance of the SLOCC Kirwan polytope from the 
origin. 

1. Introduction 

The problem of classification of pure multipartite entanglement is a recurrent subject 
in the quantum information theory. Here we analyze two possible approaches to the 
problem. The idea behind the first one is based on a perhaps trivial but very fruit- 
ful observation that entanglement has to be invariant under local unitary operations, 
i.e. unitary operations applied independently to every subsystem of a given system. 
In terms of group theory this amounts to the study of orbits of the compact group 
K = SU{N 1 ) x ... x SU(N L ) on the compact manifold M = P(H), where H = 
Hi <S> ■ ■ ■ <£> T~Ll is the Hilbert space of the whole system and P(H) is the corresponding 
complex projective spaceQ Since all such orbits are closed the quotient space M/K is 
at least a topological Hausdorf^] space and one can perform the desired classification 
by searching for if -invariant polynomials. 



For indistinguishable particles we need a slight modification of this setting as it will be explained in 
Section]?] 

2 A topological space is Hausdorff if any two points can be separated by open sets. 
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The second approach uses the fact that entanglement does not increase under invert- 
ible local transformations. Relaxing the unitarity condition produces a class of transfor- 
mations known under the name of stochastic local operations assisted by classical com- 
munication (SLOCC). The SLOCC operations are represented by G — SL{N\, C) x 
. . . x SL(Nl, C) and G = K c , i.e. the group G is the complexification of the group K. 
In the language of the group theory this means that G is a reductive group. However, 
since G is not compact its orbits in M may not be closed and hence the space M/G is 
typically not Hausdorff. Nevertheless, the celebrated theorem of Hilbert and Nagata ifll 
ensures that the ring of invariant polynomials is finitely generated and it is possible to 
distinguish between different closed G-orbits in M. The construction of quotient space 
M/G is however more elaborated. It is known under the name of geometric invariant 
theory quotient (GIT) which we discuss in Section[3] 

The space of states M = P(K) is obtained by identification of vectors from T~L 
differing by a multiplicative non-zero factor, 

H 3^ H € V(H). (1) 

Complex projective space is a Kahler manifold and hence a symplectic manifold with 
the Fubini-Study symplectic form ujfs- Although M is symplectic the restriction of 
ujfs to an orbit of K is typically not symplectic. Inspired by this observation we 
have recently introduced a discrete entanglement measure defined as the degree of 
degeneracy of uipg restricted to if -orbits [2|. In simple words we showed that more 
non-symplectic K -orbit through the state [v] is, more entangled is the state [v]. The 
key ingredient in the proof of this statement was the concept of the momentum map 
fj, : F(H) — > t*, where t* is the dual space to the Lie algebra t of the group K. In 
our setting this map encodes the information about expectation values of the local ob- 
servables, i.e. it is given by the collection of reduced one-particle density matrices. The 
degree of the degeneracy of ojps at M is the dimension of the set of if -equivalent 
states which have the same reduced one-particle density matrices. This degeneracy can 
be thus interpreted as a 'measure' of an additional information which one has to gain in 
order to distinguish state [v] from the states which are if -equivalent to [v] using only 
one particle density matrices. 

The idea of the momentum map proves to be fruitful also for characterization of 
G-orbits. One of the interesting and, as we show, useful properties is that the ordered 
spectra of reduced one-particle density matrices corresponding to the states from the 
closure of G-orbits form a convex polytope. This observation is in fact a direct corollary 



from the momentum map convexity theorems which we discuss in details in Section 3.2 
In the current paper we focus on the investigation of SLOCC entanglement classes for 
pure states. Our description is through the critical sets of the total variance of state, 
Var[i>], which was recently introduced by Klyachko as a measure of entanglement [3]. 
Using the Riemannian metric on F(H) we construct the gradient vector field of Var[w]. 
The limit sets of the flow generated by this field are critical points of Var[i>]. Making use 
of the momentum map geometry we notice that the function Var[w] is up to a constant 
the square norm of the momentum map ||/i|| 2 . The critical points of the later were 
studied in early 80's by Kirwan and Ness 00. Translating their results into language 
of entanglement we find that each critical set of Var[w] can be characterized as a unique 
if -orbit inside G-orbit on which Var[v] attains maximum. Therefore moving the state 
along the gradient path always increases the total variance associated with it. Moreover, 
the gradient of Var[w] is tangent to G-orbits, and hence its flow can be realized by local 
operations. In order to describe the structure of SLOCC classes, i.e. the structure of 
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G-orbits we identify all states which are carried by the gradient flow into one critical 
if -orbit. These might be not only all states from the G-orbit that contains the chosen 
critical K -orbit, but also those which can be taken to it by the gradient flow in the limit. 
Taking into account the physical meaning of Var[t>] it seems to be natural procedure. 

In order to illustrate the above described idea we consider several examples for both 
distinguishable and indistinguishable particles. The main difficulty met during calcu- 
lations is an effective way of finding critical sets of Var[u]. As we explain throughout 
the paper this problem can be divided into two essentially different parts. The first one 
boils down to finding all states for which reduced single-particle density matrices are 
maximally mixed. We describe how to do it in sections 4 and 5. The second problem 
concerns critical sets of Var[i;] among states whose G-orbit in the Hilbert space H 
contains the null vector. We show that it can be greatly simplified by using the convex- 
ity property of the momentum map. Moreover we also notice that the value of Var[v], 
or actually square root of it, calculated in the critical point can be understood geo- 
metrically in terms of distance of the SLOCC Kirwan polytope from the origin. This 
introduces a kind of hierarchy among different SLOCC classes (or more generally fam- 
ilies of classes). Finally, we point out that the Morse index at a critical point have a 
nice physical interpretation of number of non-local directions in which entanglement 
increases. In the course of argumentation we exhibit a connection between the norm 
square of the momentum ||^|| 2 (and a fortiori the total variance of the state Var[u]) and 
the second order Casimir invariant for the group K that allows for an identification of 
them as an expectation value of an observable and making them, at least in principle, 
measurable in experiments. 

Recently, after completion of our work, we have been informed that Walter, Do- 
ran, Gross, and Christandl ll3D are independently studying the SLOCC Kirwan poly- 
topes, which they call entanglement polytopes in the context of classification of SLOCC 
classes of states. Their work uses similar methods to ours and is very interesting. How- 
ever, we would like to note that there are also many conceptual differences between 
both approaches. 



2. Momentum map, complex projective space ¥(W) and G-orbits 

Let K be a compact connected semisimple Lie group. Let us assume that K acts on a 
symplectic manifold (M, uj) in the symplectic way, i.e. leaving the symplectic form ui 
invariant, 

$:K X M -> M, <P g = $(g, •), = lj. 

For such an action there exists a momentum map |4|, i.e. the unique map, jit : M — > I*, 
from the manifold M to the dual space {* of the Lie algebra { = Lie(K), such that 

1. The map /i is equivariant, i.e. /u(<? g (x)) = Ad*/i(x), where Ad* is coadjoint action 
(Ad*a, £) = (a, Ad s -i£) = (a, <7 _1 £#), where ( •, •) is the pairing between { and 
t*, and Ad g denotes the adjoint action of K on {, i.e. Ad s £ = g£,g _1 for g £ K and 
i- I"- 

2. For any £ € t, the fundamental vector field, 



t _$exp(t(,) 0) (2) 



is the Hamiltonian vector field of the Hamiltonian function fi^(x) = (n(x), £), i.e. 
= •)■ 



4 



Adam Sawicki, Michal Oszmaniec, Marek Kus 



2.1. Momentum map on the complex projective space. For our purposes we need the 
momentum map for the complex projective space, P(T~L). Before presenting its con- 
struction let us review symplectic properties of F(H). Let G be a complex reductive 
grourj^] i.e. it is the complexification, G — K c , of its maximal compact subgroup K. 
Let us assume that G acts linearly on W ~ C N . For simplicity of presentation we will 
assume that G is actually a subgroup of SL(H, C) and hence K can be chosen as a sub- 
group of SU(H). We denote by (-|-) the if-invariant scalar product on H. The action 
of G on H induces in a natural way an action of G, and hence of K, on the complex 
projective space P(T-L), g[v] = [gv]. In the notation we took advantage of the linearity 
of the group action on W and we will consequently use gv instead of in this case. 
For further reference let us make two remarks. 

1. Since the group SU(H) acts transitively on P(W) the tangent space T[ V }P(H) at [v] 
is spanned by the fundamental vector fields |2]) with the above described action of 
SUCH), i.e. $ exm) [v] = [e*H £ f su(H) = Lie(SU(H)), 

2. The tangent space T[ V ]F(H) is obtained by pushing forward vectors from T V W by 
the projection ([TJ. Let thus 

d 

x = — 
dt 



v(t) (3) 



be the tangent vector to the curve t H> v(t) 6 H at v(0) = v. The corresponding 
projected curve in ¥(T~L) is given as 

, v(t) , v I v(t) y 

\\v(t)\\ 

Differentiating we find the corresponding tangent vector z to P(H) at [v] 

Z = R-R<» (4) 

Finally let us remind that the projective space F(H) is a Kahler manifold. In particular 
it means it is a symplectic manifold. The symplectic form at a point [v] € F(H) is given 
by the formula |4) 

w (i>6 = ^t-tt , Ci) 52 G su(K), (5) 

where ^ and ^2 are the fundamental vector fields associated to £ x and £2 wa ® (see 
the first remark above). At the same time F(H) is a Riemannian manifold with the 
Riemannian metric 

0(?l,?2j = ——2 ■ (O) 

Both structures are mutually compatible and compatible with the complex structure, i.e. 
6(^1 , £2) = ^2) which makes P("H) a Kahler manifold. 

Since (• | ■) is if -invariant and the symplectic form cj is given in terms of it, the action 
of K on P(H) is symplectic. This in turn implies that if K is semisimple (which is the 
case at hand) then there exists the unique momentum map ji : F(H) — > t* . Now we are 
ready do describe 



In Section|3]we give a general definition of a reductive group from which the present description follows. 
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The construction of the momentum map. For any ueWwe define the map 

n v : G -> R, n v (g) = (gv\gv) = \\gv\\ 2 

The derivative of n v at the group identity, e e G, is by definition a linear functional 
dn v : q — > R and hence an element of g* . Notice however that since (• | •) is if -invariant, 
for any ^e!we have 

(dn v , = 4 (cxp(«>| exp(tC)w) = 
at t=o 

Next, since G is reductive, the Lie algebra g = t ® it and in fact dn v e it* . We define 
the map 

1 nu 

Straightforward calculations give 

(M,i0 = l { ^-,^t, (?) 

It is easy to see that for fi^ = i£) one gets cfy^ = •). One can also check 

that /j, is equivariant, i.e. 

MH) = Ad>(M), ff6 Jf. (8) 

In order to avoid confusions we clarify that the Hamilton function associated to the 
Hamiltonian vector field £, where £ e { is /i^. Equivalently, one can define the true 
momentum map p,^ = ^ which has all required properties and in addition satisfies 
dp,^ = •). In the following we will use the former convention. 

3. The function ||^|| 2 and G-orbits 

Let us concentrate on properties of the norm square of the momentum map on a set of 
its critical points. First, notice that for practical reasons it is customary to identify it* 
with it by means of the Ad^ invariant scalar product (-|-) on t. More precisely, for any 
a e it* we define a* <G it such that 

(a,0 = («*IO- (9) 
To give an explicit formula for fi* (x), where x = [v], let us choose orthonormal basis 
{^i} of it. Then 

dim K dim K dim K 

six) = e fa* = E (Ms), &>& = E kh( x )& 

k=l k=l fe=l 

After such an identification we define 

|2 



HmII 2 :P(W)-^R, 

dim K 1 dim K ✓ , , . \ \ 2 

( x) = = E = J E ( ^ ) ■ < n > 

fe=i fe=i 



Let us note that since \i is if -equivariant the function | | 2 is if-invariant, i.e. 

IMI 2 (M) = IImII 2 (M), geK. 
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Physical meaning of ||/i|| 2 . Remarbly, we can give to \\/j,\\ some "physical" meaning 
in the sense that it can be expressed as the mean value of an observable. We define after 
Klyachko [3| the variance of a state [v] 6 P(H), with respect to the symmetry group 

K c SU(H) : 



i K „ dim K 



The above expression is simply the sum of the variances of the observables & calculated 
in the state [v] . On the other hand, 

dimJf 

i=i 

is the representation of the second order Casimir operator Q, which belongs to the 
center of the universal enveloping algebra il(t) of f and therefore commutes with each 
£i . We can now rewrite Equation ( 13 1 in the following manner 

Var(H) = \{v\C 2 \v) - £ («l&l«> 2 j • (13) 

If the group X acts on P(%) irreducibly, which is always the case in our setting, then 
C2 is proportional to the identity operator and thus ( 13 1 boils down to 

/dim k \ 

Var(H)= C --— 5 E HfOI«) a =c-4.|N| 2 (M). ( 14 > 



i=l 



{ V \ C \u 1 2 

where c = v * s a M -independent constant. Notice that since ||/i|| ([«]) is K- 

invariant, the variance Var([w]) is also if -invariant and thus can serve as an entangle- 
ment measure. This argument relates the norm square of the momentum map with the 
entanglement measure considered in Q. Alternatively ||/^|| 2 ([v]) can be expressed as 
the expectation value of C2 represented on the symmetric tensor product T~L V H (see 
also J6l). Indeed, one easily checks that for 



dim if 

C 2 V = E J + ^®6) 2 , (15) 

i=l 



which is the second order Casimir operator on % V % and we have 

^ ^ dim K 

'v®v\C^\v®v)=2c + 2—- 1 V (u|&H 2 = 2c + 8|| M || 2 (M),(16) 



<^} 2 " ( v \ v f i= i 



where c is the number that appeared in ( |14) , Finally, let us notice that by the formula 
( 14 1 the critical sets of the total variance of state function are exactly the critical sets 
of~[|/i|| 2 . In the subsequent sections we show how one can use these critical sets to 
parametrize all SLOCC classes of states. 
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3.1. Critical points o/||/i|| 2 . It was noticed around quarter century ago [7,8] that crit- 
ical points of ||/i|| 2 play an important role in the classification of G = i\T c -orbits in 
F(H). In the next two paragraphs we give the description of the set of critical points of 
||/i|| 2 . Let us start with a 



Formula for d\ \fi\ \ 2 (x). From Equation ( 11 1 we have 



/dim if 



dim if 



dim if 



d\\iA\ 2 i x ) = d [ Y t4i( x )\ = 2 Y ^( x ) d ^ = 2 Y1 MfJ^M-^fe. 



i=l 



dim if 



i=l 
dim if 



2 ^(/%0*0-*£fc,-) = 2u>(-i Y Vik( x )€k,-) = 2u)(-ifi*(x),-) = 2d/i ft .( x ). 

(17) 



fc=i 



fc=i 



Notice also that since the symplectic form cj is nondegenerate rf||/i|| 2 (a;) = if and 
only if — ifx* (x) = 0. 



Two kinds of critical points of \\(J,\\ ■ Using formula (17i we can divide the set of 
critical points of | | 2 into two disjoint classes 

- The so called minimal critical points for which fi(x) = 0. Notice that if /i(x) = 0, 

or equivalently fJ,*(x) — then of course — ip,*(x) — and hence c?||/i|| 2 (a;) = 0, 
i.e. a; is a critical point. All points in the jti _1 (0) are thus critical. 

- The non-minimal critical points for which fi(x) ^ and — ipj* (x) = 0. 
Let us observe that [v] is a critical point if and only if 



H*([v])v = Xv 



(18) 



for some A G C. Indeed, using ^ for the curve t (-> exp(— ip,*([v])t)v we obtain 
x = -ifi*([v])v and from (p} 



-*M*([«]) = rr^ U*([v])v- j^(v\l**([v])v)v 



(19) 



Hence, —ifi*([v]) = implies ( 18 1 with 



(«|m*(M)«) 



(20) 



Conversely, if (fT8~| l is fulfilled then by scalar multiplying its both sides by v we find that 
A is given by (20i. Consequently, upon \l9\ we obtain —i(i*([v]) = 0. 
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The minimal critical points and the construction of the G-orbit space. . The minimal 
critical points of | | 2 play an important role in the construction of the G-orbit space in 
F(H). For simplicity we begin with the construction of the G-orbit space on the level 
of H. Notice first, that since G is not compact its orbits in T~L may not be closed and 
in general the topological quotient space H/G is not Hausdorff. Our aim is to describe 
how to overcome this difficulty and construct a "good" quotient space. We first recall 
the most basic definition of a reductive group. To this end let V and U be two complex 
vector spaces. Let <f> v : G — > GL(V) and <\F : G — > GL(U) be two representations 
of G on V and U respectively. The so-called intertwining operator is the linear map 
A : U — » V, such that for any vector u € U 

4% oAu = Aocf>V(g)u. 

We denote by U G and V G the sets of G-invariant vectors, i.e. the vectors in U such that 
4>gU = u and the vectors in V such that <f>g v — v for any g 6 G. The operator A maps 
U G into V G . To see it notice that 

V<7 G G 4> u g u = u=> (f^ o Au = Ao 4>U{g)u = Au o Au e V G 

We are now ready to give the definition of the reductive group [ 1| . 

Definition 1 An algebraic group G is reductive if and only if for every surjective inter- 
twining operator A : U —> V of G -representations, the induced map on G-invariants 
A : U G —¥ V is surjective as well. 

Going back to our setting we have the natural action of a reductive group G C SL(H) 
on H. This action induces an action (representation) of G on the ring of homogenous 
polynomials C[H] 

9-f(v) = f(g- 1 v). 
The following theorem holds for reductive group actions. 

Theorem 1 (Hilbert, Nagata [1]) If G is reductive, then the ring of G-invariant poly- 
nomials C [H], i.e. polynomials which are constant on the G-orbits is finitely generated. 

Having Theorem [T] one can ask if it is possible to distinguish between all G-orbits 
in H using invariant polynomials. The answer is unfortunately typically negative. It 
is because G-invariant polynomials are continuous and they take the same value on 
orbits for which the intersection of closures in not empty, G.V\ H G.v-i ^ 0. We 
call two orbits G.v and G.u closure equivalent if there exists a sequence of orbits 
G.v\ — G.v, G.V2, ■ • ■ , G.v n — G.u such that G.Vk fl G.vt+i ^ 0. This equiva- 
lence relation divides G-orbits into closure equivalence classes. The proper question is 
hence the following one. Do G-invariant polynomials separate the closure inequivalent 
orbits? We have the following U) 

Theorem 2 Let G.v\ and G.V2 be two orbits such that G.v\ fl G.vi = 0. Then there is 
a G-invariant polynomial which separates G.v\from G.vi- 

The immediate corollary from this theorem is 

Collorary 1 The following are true: 

1. There is enough G-invariant polynomials to separate closed G-orbits. 
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2. Closure equivalence class of orbits can contain at most one closed orbit, as if it 
contained more than one then invariant polynomials would not be constant on it. 

3. Every closure equivalence class of orbits contains exactly one closed orbit ( the orbit 
of lowest dimension). Moreover, it is contained in the closure of every orbit in the 
equivalence class. 

The last point implies that we can simplify the definition of the equivalence of two 
orbits G.v and G.u to CLv H CLu ^ 

The correct quotient construction is thus with respect to closed orbits. In the case 
of the projective space P(H) one first removes the so-called null cone N which is 
the closure equivalence class of orbits which contains £ Mas the unique closed 
G-orbit. What is left is denoted by P(H) SS and called the set of semistable points. 
Then we say that two points x\ = [vi] £ F(W) SS and x 2 = [1*2] € P('H) SS are 
equivalent, [vi] ~ [v 2 ], if and only if G.V\ n G.v 2 7^ 0. The resulting quotient space, 
denoted F(H) SS // G, is known in the literature as the categorical quotient. Notice that 
¥(H) SS //Gis parametrized by closed G-orbits. Remarkably, the minimal critical points 
of I | 2 parametrize all closed G-orbits in F(H) SS . It was proved by George Kempf and 
Linda Ness fg] that 

Theorem 3 Assume that [v] G /i _1 (0), i.e. [v] is a minimal critical point of ||/i|| 2 : 
P("H) — > R. Then the orbit G.v is closed and ^ G.v. Moreover all closed G-orbits in 
¥('H) SS are of this type. 

It is now easy to see that we have the following isomorphism of quotients 

A T 1 (0)/Ji-£*P(«)„/G, 

Moreover, we have [7] 

Theorem 4 The sets V(TL) SS // G and hence /i -1 (0) / K are projective varieties. 

By Theorem|4]one can view the categorical quotient construction in terms of fibers 
of certain surjective map between complex spaces, i.e. we define 

tt: ¥(H) ss -+P(H) ss //G, 

where each fiber of tt is a closure equivalence class of orbits containing exactly one 
closed G-orbit which in turn contains the i^-orbit from /i _1 (0). It is natural to ask 
when a fiber of it is given by a single G-orbit, i.e. when a closure equivalence class 
of orbits is simply one closed G-orbit. In order to address this question we need the 
following general theorem 

Theorem 5 Let tt : X — > Y be a surjective map between complex spaces and assume 
that each fiber of tt is connected. Let d — dim X — dim Y. The following are true 

1. There exists Zariski open dense set Y reg C Y such that for y G Y reg the fiber tt^ 1 (y) 
has dimension dim7r _1 (?;) = d. 

2. For each y £ Y \ Y reg the dimension of fiber 7r _1 (y) cannot drop down, i.e. 
dimTT^ 1 (y) > d. 

On the other hand, it is known that 

Theorem 6 The set P(H) S = {G.x : dimG.x = dimG and x E ^i _1 (0)} if exists is 
Zariski open dense subset of¥('H) ss . 
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We will call the points from F(H) S stable. Notice now that if a; € P(H) S then the 
fiber of ir : P(W) SS —> F(H) SS // G is exactly G.x. To see it let us assume on the 
contrary that there is a point y E P(W) SS such that G.y is closure equivalent to G.x. By 
Theorem [T] this means that either G.y = G.x or dim G.y > dim G.x. But the later is 
impossible since dim G.x = dim G. Therefore the generic dimension of a categorical 
quotient fiber in the presence of stable points is d = dim G. Combining this observation 
with Theorem [5] we arrive at 

Theorem 7 For a semistable but not stable point x € /i -1 (0) the fiber of categorical 
quotient n~ l (n(x)) contains more than one G-orbit. 

One would naturally want to know at this point how to find the remaining, i.e. 
non-closed G-orbits in the fibers of 7r _1 (7r(a;)) for a semistable but not stable point 
x £ /x _1 (0). The general procedure to do it is known in algebraic geometry under the 
name of the Richardson's proof of Hilbert-Mumford criterion (T3). In short words it 
boils down to construction of some particular one-parameter subgroup that carries a 
semistable but not stable point y asymptotically into X -orbit in /x _1 (0). We will not 
describe here this subtle construction. We address this problem in its full generality in 
the subsequent publication ifTUll . In Section[5] however, we give a flavor of it by showing 
that most of the SLOCC classes for four qubits introduced in [24 1 can be understood 
using this procedure. 

Finally, let us point out that it is known (see JTT)) that the set G.fi~ 1 (Q) is either 
empty or open and dense in F(H). This means that n~ l ((S)/K parametrizes almost all 
G-orbits in F(H). Moreover, all these points are in one-to-one correspondence with 
the minimal critical sets of \\fJ,\\ 2 . A natural question now is whether it is possible to 
stratify the null cone by the closure equivalence classes of orbits via critical sets of 
||/i|| 2 . In other words we want to introduce an additional refinement in the null cone iV 
stemming from the critical sets of the total variance of the state restricted to it. In the 
next paragraph we discuss the structure of these critical sets. 



G.x 2 G.X! 




Fig. 1. The idea of the categorical quotient construction, i.e. fi (0) / K = ¥(H) S s II G. 
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The Morse index at a critical point . An important information about the critical point 
of a function is encoded in its Morse index, i.e. the number of negative eigenvalues of 
the Hessian at the point. To calculate it for critical points of | | 2 let us, as previously, 
choose an orthonormal basis {£,*} of it. We already know that 

(dim K \ dim K 

J2 nU x n = 2 E Kite)***- 
i=l / i=l 

Let us assume that a; is a critical point of ||/i|| 2 . The Hessian of ||/i|| 2 (x) is the matrix 
of the second derivatives (or equivalently the corresponding quadratic form) given by 

dim K dim K 

Hess|H| 2 (:r) = 2 ^ d^ k ® dn ih + 2 ^ /^Hess/%. (21) 

fe=l k=l 

For a minimal critical point x, i.e. p,{x) — we have /i£ fc (x) — for all the formula 
pi) simplifies to 



dim/f 

Hess||/i|| 2 (x) = 2 dfi£ k (3 d[j,£ k , 

k=l 



and the Morse index, i.e. the number of negative eigenvalues of Hcss||/i| | 2 (a;) is equal 
to zero, ind(x) = 0. 

Assume now that £ is a non-minimal critical point, i.e. fi(x) ^ and — ip,* (x) = 0. 
We can always choose an orthonormal basis £j of it in a such a way that £i = n^feyn • 
This means that all /i^ fc (x) = (p(x), = (t;k\li*(x)) = for k = 2, . . . , dim If and 
Ki( x ) = (ra^'W) = The formula (E 



simplifies to 



dim if 

Hess||^|| 2 (a;) = 2 ^ d/^ fc ® d// ?fc + 2||/x||(x)Hess/X£ 1 . 
fc=i 

We have following theorem Q: 

Theorem 8 Assume that x is a critical point of | | 2 : P(%) — > R. 77zen the following 
is true 

1. \\^\\ 2 restricted to G.x attains its minimum value at x. 

2. ||ytt|| 2 restricted to G.x attains its minimum value on a unique K-orbit, which is of 
course the orbit K.x through the point x. 

3. If x G /i _1 (0) then the Morse index ind(a:) = 0. 

4. If x ^ /x _1 (0) then ind(a;) may not be zero. 

Notice that 3 and 4 have been already proved. The proof of 1 is a direct consequence 



of the convexity property of the momentum map (see Section 3.2 below). For the proof 
of 2. see Q (Theorem 7.1). Moreover, we can decompose the tangent space T X P(W) as 
the direct sum 

T X V(H) = T x G.x © (T x G.x) ±UJ . 

It is possible since T x G.x is a complex and hence a symplectic vector space and hence 
T x G.x n (T x G.x) ±ul — 0. Moreover these subspaces are also orthogonal with respect 
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to the Riemannian metric &(■, •) (|6j as they both are stable with respect to multiplication 
by i. By Theorem|8]we know that the Morse index at any critical point is determined 
by the restriction of Hess| \ 2 (x) to the subspace (T X G .x) ±ul . Notice next that for any 
vector v E (T x G.x) u we have 

dp, £k (v)is = u(£ k , v) = 0, 

and hence 

(Hess||Ai|| 2 (ar)) \ {Tj:G . x) ^ = 2||^||(x)HessM^- 

Summing up the Morse index at the critical point is given by the number of the negative 
eigenvalues of Hess/i^ where £1 = n^feyn ■ 

3.2. The convexity property of the momentum map . As the last tools for our applica- 
tions we need some facts concerning the convexity property of the momentum map. Let, 
as previously, (M, ui) be a symplectic manifold endowed with the symplectic action of 
a compact connected group K. Let us fix T C K, a maximal torus in K. As before, the 
Lie algebra of K will be denoted by { and its Cartan subalgebra (the Lie algebra of T) 
by t. Positive elements in t form the so-called positive Weyl chamber t+ |[T2l . Finally 
let p : M — > t* be the momentum map. Notice that since /i is equivariant, any K orbit 
in M is mapped onto a coadjoint orbit in J* . Moreover, every coadjoint orbit crosses t* + 
in exactly one point. We can thus define : M —> t* by <P(x) — /j,(K.x) n t+. By 
restricting the identification (|9| to t* we may alternatively treat \P(x) as an element of 
tf whenever it is more convenient. 
The following is true 

Theorem 9 The set Hf{M) is a convex polytope. 

This theorem has a long history and many people contributed into its final form. First, 
Atiyah [14] and Guillemin and Sternberg lTT5l showed that for an abelian K even more 
is true, namely the whole /i(M) is a convex polytope. Guillemin and Sternberg also 
noticed that for a non-abelian K the set n{M) is typically not convex and proved that 
!?(M) is a convex polytope when one assumes in addition that M is a Kahler manifold 
|16|. Finally, Kirwan J8J showed that Theorem [9] holds for any connected symplectic 
manifold. We will call \P(M) the Kirwan polytope. 

In this paper we are especially interested in the setting when M = ¥(%) together 
with the action of a complex reductive group G = K c . In this case the finer results are 
known. We give here one of particular importance for our purposes: 

Theorem 10 Let G.x be an orbit ofG through x € F(7i). Then 

1. The set ^(G.x) is a convex polytope SB, STB- 

2. There is an open dense set of points in ¥{%) for which tf'(C.x) = ^(P('H)) SIEl- 

3. The collection of different polytopes ^(G.x) where x ranges over M is finite IfWj . 
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3.3. Stratification of the null cone by critical sets of | | 2 . After introducing all tools 
we are ready to describe the stratification of the null cone N by critical sets of ||/x|| 2 
and provide the method for finding classes of G-orbits. Let us first notice that since 

Hamiltonian vector field corresponding to d (— | | 2 ) is given by 2i/i*(x) and P(H) is 
a Kahler manifold it means that gradient of | | 2 with respect to the Riemannian metric 

&(■, ■) is given by 2/i* (x). Next, the vector field [i*(x) is tangent to the G-orbit through 
x as /i* (x) e it C g. Thus any point x £ N is carried by the gradient flow into some 
critical A'-orbit of \\fJ,\\ 2 and in this way we obtain stratification of N. More precisely 
for each a such that a = 'P(x) with x e N, we check if 'P^ 1 (a) contains | | 2 -critical 
if -orbits. We denote by C a the set of ||/i|| 2 -critical points such that >P(K.x) = a and 
by N a all points in N for which C a is the limit set of gradient flow. Next we say that 
two points in N a are G-equivalent if and only if they are carried by gradient flow to the 
same critical A"-orbit in C a . Under this equivalence we get the desired stratification of 
N. It turns out that the set C a /K has a nice structure, namely Q 

Theorem 11 The set C a /K is a projective variety. 
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In this way the quotient of F(H) by G is decomposed into finite number of disjoint 
projective varieties 

¥(n)/G = \JC a /K 

ot 

The main one is given by Cq/K ~ p~ 1 (0)/K or equivalently ^ r ~ 1 (0) / K . The other 
are C a /K where C a is the set of all critical if-orbits in <? _1 (a). The fact that there is 
finite number of these varieties is a direct consequence of point 3 of Theorem [TU] 

An algorithm for finding classes of G-orbits. Summing up, in order to find the above 
described decomposition of F(H) one has to 

1. Find the set /j,~ 1 (0)/K. Usually it is a hard problem and requires knowledge of a 
canonical form of x € F(H) under if -action. 

2. Find sets C a /K where C a is the set of all critical if -orbits in if r_1 (o;). This happens 
to be relatively easy task as we show in the next section. 

Finally, let us note that a priori it is possible that some orbits in N do not contain 
critical points of Var. In our examples we did not observe such a situation. If they are 
however present they can be treated in the similar way as in the case of semistable points 
[7 19 1, i.e. by constructing appropriate one-parameter subgroup which carries them to 
the critical if -orbits. We will discuss this problem in full generality in the subsequent 
publication 1 10 1. 



4. The SLOCC classes for distinguishable particles 

We are now ready to apply the mathematical tools described in sections 2 and 3 to the 
characterization of SLOCC classes of states. 

Entanglement in systems of distinguishable particles. In the following we put % = 
C N <E) . . . <£) C N which is the Hilbert space of L identical but distinguishable particles. 
The group G is a direct product of L copies of SL(N, C), i.e. G = SL(N, C) xL and 
K = SU (N) xL . We have G = K c and hence the group G is reductive. Moreover, K is 
semisimple compact and connected. As mentioned in the Introduction the group G rep- 
resents the so-called SLOCC operations whereas K stands for local unitary operations. 
The actions of both G and K on H are given by 

(U u U 2 ,..., U L ).v = Ux <g> U 2 <g> . . . <g> U L v. 

The Lie algebras of G and K are direct sums of L copies of sl(N) and su(N) respec- 
tively, i.e. g = sl(N)® L and J = su(N)® L . The corresponding actions of g and f on T~L 
read 

(6, 6, ■ • • > Cl)v = (6 <8 In ® • • • <8 In + In ® 6 <8 • • • <8 In + • • • + 

where ijy is the N x N identity matrix. The momentum map is p : F(H) — > it is given 
byEl 

MM) = Oi(M) - ^In, P2([v]) - ~In,---,Pl{[v]) - -*t t In), 
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where /Oj([w]) is the i-th reduced single-particle density matrix of the state A. The 

points [v] £ F(H) for which p([v}) = are thus precisely the states for which all 
reduced single-particle density matrices pi are maximally mixed. Moreover the map 
^ : F(H) -> it+ is given by 

#([«]) = (pi(M) ~ MM) - }jIn, ■ ■ .,Pl(M) - ^I*r), 

where &([?;]) = diag (p[, . . . ,p N ), p\ > p\ +1 , and Spect(p 4 ([w])) = {p\, . . . ,p N }. 
In what follows we will need the explicit formula for pb* : 

p*([v]) = ypi([v]) - j^lN^j ®In® ...®In + In® [P2([v]) - jj 1 ^) ®In®---®In- 
+ . . . + I N <g) . . . (gs I N ® \ p2([v]) - ^In) ■ (22) 

The image of & is by Theorem [9] a convex polytope. For the system of L qubits, i.e. 
when N = 2 this polytope is given by particularly simple inequalities, namely we have 
the following theorem l20l 

Theorem 12 For a L-qubit system all constraints on the 1-qubit reduced density ma- 
trices of a pure state are given by polygonal inequalities 

for the minimal eigenvalues pi of 

Let us now give some details about finding SLOCC classes for identical distinguish- 
able particles. First we have the following theorem which gives canonical form of the 
state up to local unitaries l2~D : 

Theorem 13 For a general state & = J2 \h ■ ■ - ib) in U = <C N ® . . . ® C N 

there exist local unitaries E/j such that all the following entries in the state Of' = JJ\ ® 
. . . ® Vifb are set to zero: 

vi < j < n vfc > j q 7 ,.. ifc = c' :jj „_ jkj ... c' kJ „, j:j = o. 

Moreover all entries C' NNiN N , where i < N can be made real and positive. If the 
number of parties exceeds 2, then the normal form is typically not unique up to permu- 
tations, but there exist a discrete number of different normal forms with the aforemen- 
tioned property. 

The algorithm from Section [373] thus reads 



1. Use Theorem[l3]to find ^ _1 (0) /K. This gives family of SLOCC classes in F(U) SS . 
We show in next section how to do it explicitly for two, three and four qubits. 

2. For SLOCC classes in the null cone N we have the following strategy. Since crit- 
ical sets of \\p,\\ 2 are if -invariant it is enough to check when p,* ([v]).v = Xv for 
//*([«]) G !f - (P('H)) \ {0}. In other words we have to go over the Kirwan polytope 
and verify for which states v we have a*v — Xv, a* £ !f r (P(H)) \ {0}. For many- 
qubit systems this polytope is described by polygonal inequalities (see Theorem 12 1 
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which are explicitly known. Moreover, we can treat a* as an N L x N L diagonal ma- 
trix acting on H. If the diagonal elements of a* are nondegenerate the corresponding 
eigenvectors are separable states. In order to find other nontrivial SLOCC classes 
we are interested in the situation when the spectrum of a* is degenerate. Remark- 
ably, the dimensions of the degenerate eigenspaces are typically relatively small. 
For consistency, for each of such eigenspaces we choose these vectors v for which 
a = fi ([v]). In this way we arrive at sets X a = n C a , where a — )P([v]). In 

order to obtain set C a /K we still need to get rid of the remaining K freedom. We do 
it by constructing the space C a /K = X a /K a where K a is the isotropy subgroup 
of a with respect to the adjoint action of K. In the next section we show how to do 
it for two distinguishable particles as well as for three and four qubits. 



Before we give some details about application of the algorithm described in Section pJ] 
let us describe additional ideas connecting critical points of ||/i|| 2 and the convexity 
property of the momentum map with entanglement. 

1. To any point [v] E F(H) we can associate a positive number in the following 
way. For [v] £ P(H) consider the convex polytope fy(G\v\). We put d([v]) = 
dist(0, \P(G.[v])), i.e. d([v]) is the distance from to the polytope W(G.[v]). Since 
d([v]) = infygG.HllMj/)!! we nave ^(M) — f° r M e M _1 (0)- Moreover, since 
€ \P(G.x), this distance is also zero for any x G G.^ _1 (0). In other words the 
states which can be transformed by SLOCC operations to a state for which all one 
particle density matrices are maximally mixed are characterized by d([v]) = 0. 
Hence d([v]) ^ if and only if [v] is in the null cone N. 

2. Using the results described in Section [3] we get a nice characterization of points 
x G G.[v] for which | (x)| | = By Theorem|8]we know that these are points 
from the critical if-orbits of ||/i|| 2 . Moreover by definition of the total variance of 
a state Var[u], the states belonging to this critical A"-orbit are the most entangled 
representatives of the appropriate SLOCC class. 

3. The Morse index at any point [v] belonging to the critical if-orbit equals to the num- 
ber of independent directions representing the (non-SLOCC) directions in which the 
variance Var[w], or equivalently entanglement, increases. Notice that by Theorem[8] 
for states in /i _1 (0) this index is always 0, i.e. even non-SLOCC operations can not 
increase the total variance of states belonging to /i _1 (0). On the other hand for a 
separable state [v] the Morse index is maximal and equal to dimP('H) — dim G.[v]. 
It is simply because any non-SLOCC operation increases entanglement of separable 
states. 

Summing up, both d([v]) and the Morse index have an interesting geometric and physi- 
cal interpretations and can be considered as entanglement measures. Notice finally that 



by equation ( 14 1 which relates Var[u] and \ \fi\\ ([v]) we can rephrase d([v]) in terms of 



square root of the total variance 

d([v]) = \ sup„Var([u]), 

where the maximum is taken over the states [u] belonging to the SLOCC class of state 
[v]. Finally, let us clarify that by the above dicribed procedure we divide the set of 
pure states F(H) into three disjoint classes. The first one contains states which can be 
converted by G to states characterized by maximally mixed reduced density matrices. 
The second one posses the same property only asymptotically and the third one does 
not have it. Moreover, the first class, if exists, contains almost all states from F(H). 
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5. The main examples for distinguishable particles 

We will illustrate the ideas described in Section|4]by calculations for bipartite states and 
three and four qubits. 



5.1. Bipartite case. Let us consider two identical but distinguishable particles, i.e. the 
Hilbert space of the system is H = C N ® C . The SLOCC operations are elements 
of the group G = SL(N,C) x SL(N,C) and local unitary operations of the group 
K = SU(N) x SU(N). It is well know that up to a iVT-action any state v G % can be 
written in the so-called Schmidt form, i.e. 

N N 

v = y^a»|i) (g> \i) = y^^ai\i,i}, (23) 

i=l i=l 

where {li)}^! is an orthonormal basis in H, a; > a nd y^ =i af — 1. Our first task is 
to determine the critical points of ||/i|| 2 . From Section 3.1 we know that [v] € P('H) is 
a critical point of | |/.t| | 2 if and only if 

M*([tij) = IM*([v])v = Xv, 
for some A 6 C, i.e. when v is an eigenvector of 

Pi([v])-^In) ®I n + I n ® \ P2{[v}) - ^In)^) v 

where /?i([t>]) are the reduced one-particle density matrices. Notice that since ||/i|| 2 is 
^-invariant it is enough to consider states v in the Schmidt form, which we assume in 
what follows. For the state d23l we have 



and 



pi(M) = p 2 (H) = diag (a\, a 2 ., . . . , a 2 N ) , 
N ( 1 \ 



To verify when 



N / , s N 



I 1 \ 

E K- ^ H™> = A E a ^' i >' (24) 

i=l ^ ' i=l 

we consider two cases 

1. All Schmidt coefficients ai, ■ ■ ■ , ajy ^ 0. In this case the condition (24i implies that 
a 2 -l=A, foralKe {!,..., N}. 



Since a, > and Yli=i a 2 = 1 it is possible if and only if A = 0, i.e. a,; = for 
all i 6 {1, ... , N}. 
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2. The coefficients Pk+x, ■ ■ ■ ,Pn — 0, where k > 1. In this case the condition (24i 
imposes 

af - ^ = A, for alii € {1, ... , k}, 
which is possible if and only if and only if A = ^jr, i.e. pi — 4= for all i e 

{i /••!• 

Summing up, 

Fact 1 The critical points of \ \ \i \ \ 2 are K -orbits through the maximally entangled states 
of rank 1, . . . , N, i.e. states of the form 

1 k 
Vk = -7=y2\i,i), 

^ k i=l 

where fee {1, . . . , N}. By the formula ( |i4| ) the states belonging to each orbit K.Vk are 
those with the maximal variance Var([uJ) among all states in G.Vk, i.e. in the SLOCC 
class of state Vk ■ 

Once we know the critical sets of | | 2 we can easily compute the distance of each 
SLOCC class polytope ^(C^]) from the origin, i.e. the numbers d([vk})- 

Fact 2 The distance of the SLOCC class polytope W(G.[vk]) from the origin is given 
by 

d([v k }) = Jtt (pi(N) - ^l) + tr (p2([v k ]) -j^I 

__ y/2 (k(N - k) 2 + k 2 (N- kj) 
~ Nk 

where fee {1, . . . , N}. 

The next task is to calculate the Morse index for each critical set. Let us recall that 
by Section 3 it boils down to calculation of Hess^,^] restricted to (Tr Vk iG.[vk}) 
First we notice that 

(T[ Vk ]G.[v k ]) W = Span{|m,n), i\m,n) : m,n € {k+ 1, .. . , N}} . 

It is thus enough to consider a perturbed state of the form v — Vk + Ylij=k+i 
and the function 

(v\fi*([v k ])\v) (v\pi{[vk]) ® I + I ® p([vk])\v) 

/V(M) = TT~\ = TT~\ ' (25) 

(v\v) (v\v) 



Making use of the assumption that a-k+i, . . . , ajv = for v k we can rewrite (25 i in the 
form 

„ , L .n _ (^IMM)K) 2 E^Lfc+i \°ij\ 2 n „ 
(v\v) TV (v\v) 

Hence any perturbation of in a direction from (T[ Vk ]G.[vk\) ±UJ decreases its value. 
Finally thus, 
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Fact 3 The Morse index at critical set K. [vk] is given by 

md([v k ]) = 2(7V-fc) 2 , 

where k £ {1, . . . , N}. 



5.2. Three qubits. The simplest non-bipartite case involves three-qubits. The Hilbert 
space of the system is T~L = C 2 <g>C 2 <g>C 2 . The SLOCC operations are given by the group 
G = SX(2,C) x3 and the local unitary operations by the group K = SU(2) x3 . We 
know that in order to find SLOCC classes we have to find eigenvectors of the equation 
( p~8] > which for three qubit reads 

Pl - ^\ ® I ® I + I ® (f>2 ~ ^ ® / + I ® I ® - ^ ^ U = A«. 

Following the algorithm given in Section [4] we have to find all states for which all 
matrices pt are maximally mixed. To this end we use the canonical form of the state 
which by Theorem[l3]is given by 

v = p|011) + <?|101> + rjllO) + s|lll) + z|000), (27) 

where z € C and other coefficients are real and positive. The reduced one-qubit density 
matrices of state v in the form (|27]i can be easily calculated: 



s 2 



A 


2 + 


P 2 




ps 




ps 




q 2 - 


h r 2 


z 


2 + 


q 2 




qs 




qs 




P 2 - 


hr 2 


* 


2 + 


r 2 




rs 




rs 




p 2_ 


hq 2 



pi(M) = 



Ps(M) = 

We require that all diagonal elements of pi are the same and the off-diagonal vanish. It 
leads to the equations 

2 2 2 I 1 2 2 n 

p = g = r = |z| — s , ps = gs = rs = U, 

having two solutions. The first one, s = and p 2 = q 2 = r 2 = \z\ 2 , corresponds to the 
state 

vi = | (|011) + |101) + |110) + |000>) , 
and the second one is s ^ and p = q = r = 0, i.e. 

„ 2 = -L(|ooo) + |iii», 

where in both cases we have already removed the freedom of phase stemming from z. 
Notice, however, that the states v\ and V2 are local unitary equivalent] i.e. 

U®U<E)U.v 2 = vi, 



4 For yet another justification of local unitary equivalence see 1291 . where it is shown that K.V2 is La- 
grangian. 
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where 

In this way we find that the SLOCC class corresponding to /j,^ 1 (0)/K can be repre- 
sented by a single state vqhz — ^75 (1 000) + | 111)) which is the well-known Greenberger- 
Horne-Zeilinger (GHZ) state. What is left is to find remaining SLOCC classes in the 
null cone N. To this end it is enough to consider diagonal pi and we assume that 

1 - Pi 
Pi 

where < pi < \. This, of course, means that 

K-^-y % Pl -i) = {o-x l 

where \ > Xi > 0, i.e. ^*([t>]) belongs to the positive Weyl chamber. We can always 
write /i* ([v]) as the 8x8 diagonal matrix in the basis {| ijk)} 



/**([«]) = diag(/i, f 2 , f 3 , / 4 ,-/4,-/ 3 ,-/ 2 ,-/i), (28) 



where 



/1 — Ai + A2 + A3, ji — X\ + A2 — A3, fs — Ai — A2 + A3, 

h = Ai - A 2 - A 3 = -/1 + h + h- (29) 

We are interested in the structure of spectrum of /x*([i>]) on the Kirwan poly tope which 



by Theorem 12 is given by the inequalities 



- > -Ai + A 2 + A 3 , - > X 1 - A 2 + A 3 , - > A! + A 2 - A 3 , (30) 

where 1 > A^ > 0. In order to find the SLOCC classes in the null cone N we assume 
that fi*([v]) and analyze eigenspaces of n*([v}) on the Kirwan polytope. The 



following fact is a straightforward consequence of d28b and (29 1 



Fact 4 Assume that 7^ 0, then /!*([«]) has no eigenspaces of dimension 5, 6, 7 

and 8. 

Therefore we have to consider only a situation when the spectrum is non-degenerate 
or 2, 3, 4-fold degenerate. Notice, that if the spectrum of fi* ([v]) is multiplicity free then 
eigenvectors are separable vectors \%jk). However only the vector |000) gives /i(|000)) 
in the positive Weyl chamber. Hence |000) represents the first SLOCC class contained 
in the null cone. Assume now that two diagonal entries of ( |28] i are the same. There 
are a priori (Z\ = 28 possibilities to consider. However all of them lead to states for 
which either reduced density matrices are not diagonal or for which degeneracies in the 
spectrum of /j,*([w]) are of higher order. Let us now consider the case when spectrum 
of has a 3-fold degenerate eigenvalue. Notice, however, that the equality of 

any three among four /^-s implies the equality of all four f. This means that 3-fold 
degenerate eigenvalues come from choosing two /^-s and one — /fe. Going over these 
possibilities we find that the only nontrivial one is given by /2 = /; = —fi, i.e. 
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Ai = A 2 = A 3 = a and H*{[v\) = {3a, a, a, —a, a, —a, —a, —3a}. The general 
state from the eigenspace corresponding to the eigenvalue a is hence 

v = z 1 |001)+z 2 |010)+z 3 |100), 

the matrices pi — \l for [v] read 

N 2 + M 2 -| o \ /| Zl |2 + N 2_i o 

o NI 2 -|J'l o |z 2 | 2 -| 

^l 2 + N 2 -| o 
o 

The solution satisfying Ai = A 2 = A3 is \zi\ 2 = \z 2 \ 2 — 1 2:3 1 2 , so 

«=-^(|001> + |010> + |100>). 

Notice also that the state corresponding to the — a-eigenspace is LU equivalent to v and 
is not mapped to the Kirwan polytope. Finally for 4-fold degenerate eigenvalues we 
have the following possibilities 

1- /1 = fa = h = fa, in this case 

M*(M) = {Ai, Ai, Ai, Ai, -Ai, -Ai, -Ai, -Ai}, 

i.e. A 2 = A 3 = 0. It is easy to see that this corresponds to v = |0) ® ( 1 00) + |11)). 
At the same time we get that the state corresponding to the eigenvalue — Ai is |1) <S> 
(| 00) + 1 1 1) ) which is LU equivalent to v. 

2. /1 = / 2 = -/ 3 = -fi, in this case 

M*(H) = ( A 2, A 2 , -A 2 , -A 2 , A 2 , A 2 , -A 2 , -A 2 } 

i.e. Ai = A3 = which corresponds to v = 1 000) + 1 101) . The state corresponding 
to the eigenvalue — A 2 is |010) + 1 1 1 1) which is LU equivalent to v. 

3. /1 = -/ 2 = h = - fi, in this case 

M*(H) = ( A 3, -A3, A 3 , -A 3 , A 3 , -A 3 , A 3 , -A 3 } 

i.e. Ai = A 2 = 0. This corresponds to v = |000) + 1 110) . The state corresponding 
to the eigenvalue — A 3 is 1 001) + 1 1 1 1) is LU equivalent to v. 

4. fi = fa = 0, i.e. A 2 = — Ai and A 3 = 0. Now, 

fi*([v}) - {0, 0, 2Ai, 2Ai, -2Ai, -2Ai, 0, 0}. 

This, however, by positivity of Aj implies that also Ai = 0. Similarly we exclude 
cases when other pairs of fas are zero. 
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Summing up we obtained exactly six SLOCC classes of states which are given by the 
G-orbits through ll22ll 

vghz = ^ (|000) + ]1H>) , 
v w = (|10Q) + |010) + |001» , 

VO 

v Bl = (|000) + |011)) ,vb2 = A= (1000) + |101>) , 

V B3 - 4= (1000) + |110)) , v SE p = |000). 
v2 



Fact 5 77ie critical points of \\fi\ \ for three qubits are K-orbits through states Uj, By 
formula {14) states belonging to each orbit K.v^ are those with the maximal variance 
Var([u]) among all states in G.Vk, i.e. in the SLOCC class of the state Vk ■ 



We can now easily compute the distance of each SLOCC class polytope \P(G.[vk] 
from the origin, i.e. the number d([vk\). 



Fact 6 The distance of the SLOCC class polytope W(G.[vk]) from the origin is given 
by 



ro 



d([vk}) 



\ 



1 



i—l 



for v GHZ , 
for V\y, 

for VBl, VB2, V B 3 

for vsep- 



To calculate the Morse index for each critical set K.Vf. we have to find Hcss/i M * re- 
stricted to (T Vk G.{vk]) ±u ' ■ Let us first notice that (T vghz G\vghz\)' Lui — as the 
G-orbit through the state vqhz is dense in F(H) (see for example |23|). This means 
that ind(vGHz) = 0. In case of v\y one has 



(T Vw G.[v w ]) ±u = Span{|lll), i|lll)} . 

We can hence consider a perturbed state in the form v = v\y 
Simple calculations give 



/v(M) 



(vw\n*{[vw])\vw) 



(el 



ei|lll)+*e a |lll) 



2(1 



It is now clear that md(vw) = 2. Similar calculations show that for bi-separable states 
the Morse index is md(vB k ) = 6 and for the separable state ind(vgEp) = 8. 



5.3. Four qubits. In E4ll it was shown that there are nine families of inequivalent 
SLOCC classes of four qubits. In this subsection we show how to calculate the main 
family, i.e. fi^ 1 (0)/K and prove that among other there are five which are closure 
equivalent to it. 
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The main family ofSLOCC classes, i.e. [i 1 {{))/K. Applying Theorem 1 3 to four-qubit 
states one gets the normal form: 

v = zi|0000) + z 2 |0011> + z 3 |0101> + z 4 |0110) + a 5 |0111) + z 6 |1001) + z 7 |1010) + 
+a 8 |1011) + 2i9 1 1100) + aio|H01) + a n |1110) + oi 2 |llll), 

where Zi are complex and Oj are real positive. The one-qubit reduced density matrices 
are given by 



Pi 



\zi\ 2 + 1^2 1 2 + \z 3 \ 2 + \z 4 \ 2 + o§ z 2 a$ + z 3 a w + z 4 on + a 5 a 12 

z 2 a$ + z 3 a w + ZiOu + 05012 \z e \ 2 + \z 7 \ 2 + o| + \z g \ 2 + a 2 + a? x + . 

\z\\ 2 + \z 2 \ 2 + \z 6 \ 2 + \z 7 \ 2 + a| z 2 a 5 + z e a 10 + z r a n + a%a 12 



12 



/,J • z 2 a 5 + z 6 a w + z 7 an + a s ai 2 \z 3 \ 2 + \z 4 \ 2 + a\ + \z g \ 2 t u w t u n -r u, 12 

zi\ 2 + \z 3 \ 2 + \z e \ 2 + \z 9 \ 2 + af z 3 a b + z 6 a 8 + z 9 a n + a 10 a 12 

' z 3 a 5 + z 6 a s + zgciu + a w ai 2 \z 2 \ 2 + \z 4 \ 2 + a 2 . + \z 7 \ 2 + a| + ofj + a 2 2 

, z\\ 2 + \z 4 \ 2 + \z 7 \ 2 + \z g \ 2 + a\ x z 4 a 5 + z 7 a 8 + z 9 a w + anai 2 

' z 4 a 5 + z 7 a 8 + z 9 a 10 + a n a 12 \z 2 \ 2 + \z 3 \ 2 + a\ + \z e \ 2 + aj + a 2 w + aj 2 

We require all diagonal elements of p$ to be equal to each other and the off -diagonal to 
be zero. The equations for the diagonal elements give: 

a 2 ~a 2 = \z 6 \ 2 + \z 7 \ 2 ~\z 3 \ 2 -\z 4 \ 2 
al-a 2 w = \z 3 \ 2 + \z g \ 2 -\z 2 \ 2 -\z 7 \ 2 
a% + a 2 1Q = \z 4 \ 2 + \z 4 \ 2 - \z 6 \ 2 - a\ 2 

\z 4 \ 2 + \z 7 \ 2 = \z 3 \ 2 + \z 6 \ 2 (31) 

2 _ 2 
"10 ~~ a n 

The equations for the off-diagonal elements can be written as a matrix equation Av = 0, 
where 

(a 12 z 2 z 3 z 4 
z 2 a 12 z 6 z 7 

Z 3 Z 6 Oi2 Zg 
\ Z 4 Z 7 Zg Oi2 



.4 





(a 5 \ 




as 


)••- 




aio 1 







Notice that parameters in the matrix A and the vector v are disjoint and A is a symmet- 
ric complex matrix. The solution aio, a.11 , 05, as = is always a legitimate one. The 
other possibility would be detA = and v € Ker(A). It can be checked by straightfor- 
ward although tedious calculations that this possibility does not give any new solutions. 



Hence the only solution is aio, a.11 , as, as = and the equations (31 1 reduce to 



0=\z 6 \ 2 + \z 7 \ 2 -\z 3 \ 2 ~\z 4 \ 2 

o = N 2 + M 2 -N 2 -M 2 

0= \ Zl \ 2 + \z 4 \ 2 -\z 6 \ 2 -a 2 12 

0=\z 3 \ 2 + \z 6 \ 2 -\z 4 \ 2 -\z 7 \ 2 (32) 



24 



Adam Sawicki, Michal Oszmaniec, Marek Kus 



These are four equations for 8 unknowns. Writing them in the matrix form we easily 
see that solutions are given by 

1*3 1 = M, \z&\ = \ze\, \z 2 \ = \z 9 \, \zi\ = a 12 

Hence, up to normalization the states from /j, _1 (0)/i^ are given by 

v = aiflOOOO) + |1111>) + a 2 (|Q011) + |1100)) + 
a 3 (|0101) + |1010» + a 4 (|0110) + |1001)). 

It is also clear that on the projective level the reduced space /j,~ 1 (0)/K is a complex 
projective space. Summing up 

/i- 1 (0)/X = P(Ho), 

where 

Ho = Span c {ui, v 2 , v 3 , v 4 }, 

and 

di = |0000) + 1 1111) , v 2 = |0011) + |1100), 
v 3 = |0101) + 1 1010), Vi = |0110) + |1001). 

This is the class denoted in E4l by G aoc d- Notice also that when all coefficients a t ^ 
then the state is stable, i.e. dimG.[i;] = dimG. We will now show that families L a b C2 , 
La 2 b2> L a b 3 , L a4 , L a2 o ml from E?l are closure equivalent to G a bcd- The general strat- 
egy for each family is to divide the state into linear combination of the state belonging 
to G a t,cd which we denote by \v) and the rest which we denote by \w). In all cases the 
stabilizer of \v) in G contains some subgroup P(\v)) of the complex torus T c . Acting 
with the certain one-parameter subgroup P a C P(\v)) on the whole state asymptoti- 
cally we can get rid of \w). 

Family L abC2 . 

\L abC2 ) = \v) + \w), \w) = |0110), 
\v) = ^ (10000) + |1111)) + ^ (|0011) + |1100)) + c (|0101) + |1010)) • 

p(i»» - {(- e t) • a e %). a) • a e °o i« ~i>~,-*} 

Choosing the one-parameter subgroup P a C P (\v)) 

P « = {( e°«) (V e« ) ® (V e«) ( e-)^ 61 }' 
we obtain 

P a \v) = \v), Km P Q |w) = Km e - 4a \w) =0. 
And hence lim^oo F Q |L afcC2 ) e Gabcd- 
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Family L a2b2 . 

\L a2 b 2 ) = \v) + \w), \w) = |0011> + |1100) . 
\v) = a (|0000) + |1111>) + b (|0101) + |1010>) 

e- e-0) e-r e"« 



l"» = {(o 



|a = -7,/3 



By choosing the one-parameter subgroup P Q C P (|i>)) 

P « = U l) ® ("o e-0 C J l) (V £ ) ' a G 
we get 

P a | v ) = lim P a \w) = lim c~ 2a |w) =0. 

a— >oo a— >-oo 

Consequently, lim^oo P Q |i a2 6 2 ) € G a6c(J . 
Family L a03 . 

\Lab 3 ) = \V) + \W), 

\v) = a (10000) + |1111)) + ^ (|0101) + |1010)) + ^ (jOHO) + |1001) ) , 
\w) = (|0001) + |0010) + |0111) + |1011>) , 

'<!•» = { ( A ) • k A ) « ( 'J A ) 8 (o e"') l a = ^ = -7 - 



Acting with the one-parameter subgroup P a C P (|v)) 

\{ - "J (,0 < " j { e«J I ,» I 
we produce 

P Q |u) = \v), lim P Q |w) = lim e~ 2a \w) = 0. 
Thus, again lim^oo P a \L ah3 ) e G abed- 
Family L a4 . 

\L ai ) = \v) + \w) 
\v) = a (|0000) + |0101) + 1 1010) + |1111) ) , 
\w) = i|0001) + |0110) - i|1011). 

p<!»» -{(* X) « ( e ; ° s ) « (- A) « ( o «!. ) i« — -'} 

Here we choose the one-parameter subgroup P Q C P (|u)) 

3 2a \ / e a o \ / e -2a Q \ / g-a ( 

e- 2 « e-« e 2 « e 



to obtain 

and conclude that limbec P a \L a4 ) e G abcd 



P a \v) = \v), lim P a \w) = lim e~ 2 "|w)=0. 

a. — ^-oo a— >oo 
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Family L a2 o mI . 

\La 2 mi ) = \V) + \W), 

\v) = a (|0000) + |1111)), 
\w) = |0011) + |0101) + |0110) . 

By choosing the one-parameter subgroup P a c P (\v)) 

Pa = {(^ e -3aJ ® ('o" e t) ® ('a" ^0) ® ('a" e « ) 

p a |v) = Km P Q |w) = Km c- 4a |u;) =0. 

a— >oo a— too 

And hence lim Q ^oo P a |L a2 o 3ei ) € G abed- 



6. Indistinguishable particles 

6.1. Moment map for bosonic and fermionic particles. Let K = SU(N) be irre- 
ducibly represented in the space of L bosons (Ji = Sym (C^)) or L fermions {% = 
A (C^) ). In the latter case we demand that N > L. The whole formalism concerning 

the structure and the relevance of critical points of ||/i|| 2 is valid also in this case. Only 
differences one should keep in mind is that the momentum map p, : P (H) — > it and the 
corresponding mapping p* : P (Ti) — > iu (H) are given by different expressions than 
for distinguishable particles [25 1 . The momentum map is given (up to an irrelevant 
multiplicative constant) by 

V(M)= P[v]~ j^In , (33) 
where is the reduced one particle density matrix, 

P[v] = tn (\v)(v\) . (34) 

In the above expression trj denotes that the partial trace is taken over all L — 1 parts 
of the full tensor product except for the first subsystem. It can be proven that for both 
bosons and fermions the final result does not depend upon the choice of the subsystem 
which is omitted in the partial trace. Just like for the case of distinguishable particle we 
study the map \P : P (Ti) — > it + given by: 

*([«]) =PM ~ jjlN, (35) 

where pr„i is the diagonalized reduced density matrix with the ordered spectra pt v i = 
diag (Ai, A 2 . . . , X N ), Xi > A 2 > . . . > X N . 
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6.2. Critical sets of\\ii\\ 2 . Analogously to the case of distinguishable particles, in order 
to find critical sets of Var[w] (or equivalently of ||^|| 2 ), we can consider the (unique) 
intersection of the adjoint orbit trough the image of the moment map /j, ([i>]) with the 
positive Weyl chamber it + . On the practical level this amounts to consideration of the 
diagonalized ae!f(P (H)) and the corresponding operator a*, 

a* = a ® ijv <S> • • • <S> In + other permutations , (36) 

where in each of L terms there are L — 1 identity operators and one operator a. In the 
next step one has to check which eigenstates v e % of a* indeed give a under the 
momentum map: a = fj, ([v]). The index of the critical point identified in this way 

- equals when [v] is minimal ([v] € (0)), 

- is the same as the index of the Hessian of AV(M) 

([«,]) = H ^gj )w) treated as a 

function of [w]) reduced to T [v] G.[v] ±u) C T [v] P {%), where G = K c = SL(N, C) 
acts by the diagonal action on U = Sym L (C^) or U = (C N ). 



6.3. Computation of the reduced density matrix. Let us give an explicit expression 
for the reduced one particle density matrix for indistinguishable particles. In order to 
perform this computation we use the standard bases of Sym L (C^) (for bosons) and 

f\ L (C N ) (for fermions). Let {|i)} i=1 d be an orthonormal basis of C N . 
Bosonic case - the orthonormal basis is given by the vectors: 

I 7 ) = Jf IMi))®IV(2))®---®Imjv)). (37) 

aen(L) 

where II{L) denotes the set of permutations of {1, . . . ,L}, Af is the normalization 
constant, I : {1,...,L} — > {1, . . . , L} is an arbitrary map defined by i(k) = ik ■ Two 
maps /, I' with the same values and multiplicities of occurrence of each of these values 
give rise to the same state |7). 

Fermionic case - the orthonormal basis is given by the vectors: 

I 7 ) = Jf s 9n(cr)\ia(i)) ® |V(2)> O • • • ® |V(jv)), (38) 

aen(L) 

where N is the normalization constant and I : {1, . . . , N} — > {1, . . . , d} is an injective 
map defined by i(k) = ik- 

Take an arbitrary pure state of a bosonic or fermionic system: 

v = J2 a i\ I ), ( 39 ) 

lex 

where aj e C, 2~2iex l fl/ l 2 = ^ an ^ ^ s tne appropriate set (depending whether we 
consider symmetric or antisymmetric states) of generalized indices. Under the above 
introduced notation the expression for the reduced density matrix reads as: 

P[v] = EE ojojl JX J\ )=YY a /«./ tr i(IV( W )> (\I)(J\) (40) 
\iexjex ) iexjex 

lex j~i 
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where symbol "~" in the last expression denotes the fact that summation is over J such 
that the partial trace tr x (\I)(J\) does not vanish. It is easy to check that tr! (\I)(J\) will 
be when multi-indices / and J differ in more then two labels. More detailed analysis 
will be provided on specific examples. 



6.4. Two bosonic particles. Consider V. = Sym 2 (C N ) with groups K = SU (N) and 
G = SL(N, C) acting on it in a natural way. Any normalized v G T~L can written in the 
form E6): 

v = ai (|1> ® |1» + a 2 (|2) ® |2» + . . . + a d (\N) ® |JV» , (42) 

where Y^i=i | a.^ | 2 = 1, and d is a suitably chosen orthogonal basis of C d . 

The reduced density matrix is now easily computable, 

N 

p H =X>| 2 |*Xi|. (43) 



According to < fl8) >, in order to find critical points of ||/i|| we look at solutions of the 
equation: 

fi* ([v])v = Xv, (44) 

where: 

^* (H) = \ P[v] - ^n^J ®In + In® \ f>[v] - j- 
Comparing the terms we get, 



-jjlN). (45) 



^(l^l'-^j =Aoi, * = 1 > ...,JV, (46) 

for some real A. Taking into account the requirement of ordering of the spectra of 
and the normalization we get N families of (nonequivalent) critical states: 

1 k 

v k = -=Y<\ i )®\^i k = l,...,N. (47) 

In order to compute the Morse index of the critical point [vk] we notice that : 

(T[ 1Jfc ]G.[wfc]) u = Span {|m, n), i\m, n) : m, n E {k + 1, . . . , A^}, m 7^ n} . (48) 

It is important to note that the states |m, n) appearing in the above expression are sym- 
metric states (in the sense of (|37j) and therefore dim ((T [vk] G.[v k ]) ±u> \ = (N - k) (N- 

k + 1). Analogously, like in Equation |26|, all directions from (Ti Vk ]G.[vk\) u corre- 
spond to the decrease of \\n\\ 2 and therefore: 

ind([v k ]) = (N-k)(N-k + l). (49) 

Notice that the index of [v k \ differs from the case of distinguishable particles even 
though the form of the state [vk] € P (H) that encodes the k th the critical orbit of K is 
the same in both cases. 
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6.5. Two fermionic particles. Now H = f\ 2 (C N ) with K = SU{N) and G = 

SL(N, C) acting on it. Computation of critical points of \\fi\\ 2 goes essentially along 
the same lines as in the case of bosons. It uses the following decomposition of v l27l : 

v = a, (|1) A |2)) + a 2 (|3) A |4)) + . . . + j (|2 [f J - 1) A |2 [f J - 1)) , (50) 
I — I 

where Ei=i \ a i\ 2 = •"■> an( ^ {K)}j=i ;v ^ s a suitably chosen orthogonal basis of C N . 
By repeating the reasoning from the previous section we get 1^1 families of nonequiv- 
alent critical sets parametrized by the states 

1 fc 

v k = -^5^|2t-l)A|2i), fc = l,..., 

Note that (0) is empty when is not even. The index of a critical point is computed 
in the same fashion as for two distinguishable particles and two bosons. One checks 
that: 

(T [vk] G.[v k }) ±UJ = Span {\m) A \n), i (|m) A \n)) : m, n G {2k + 1, . . . , N}} . 

(52) 

It follows that dim (T[ Vk ]G.[v k ]) ±U1 = (N - 2k) (N - 2k - 1). It turns out that each 
direction from (T[ t)fc ]G.[ufe])" L corresponds to the decrease of \\fJ.\\ 2 and thus: 

ind(H) = (N-2k)(N-2k- I). (53) 



(51) 



6.6. Arbitrary number of two-state bosons.. For L two-state bosons or, equivalently, 
symmetric states of L qubits we have % — Sym L (C 2 ) with K = SU (2) and G = 
SL (2, C). We take a basis of C 2 consisting of two orthonormal vectors 1 1) and |0). An 
arbitrary state v 6 H can be conveniently written as a linear combination of the basis 
elements (called the Dicke states) of Sym L (C 2 ), 



v = ^2a k \k, L) 



fe=0 



)) ® |V(2)) 



(54) 

|V(JV)}, k = 0,...,L, (55) 



where X)fc=o \ a k\ 2 = 1- The corresponding reduced density matrix can be explicitly 
computed in the basis { 1 1), |0)}: 



P\v] 



1 

L 



T,k=o k \ak\ 2 Efc=o 



E 



W— 1 afcOfc+i 



\/(fc)(fc + l) 

fc)|a fe | 2 



(56) 
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In order to find fi 1 (0) (that is minimal critical points of ||^|| 2 ) one has to impose 
condition p^j = |l, which amounts to solving equations 

L 



I>ki 2 = 

k=0 
L-l 



j2 

fc =o \J ( L k ) ( k +l) 



L 

2 ' 

= 0. 



(57) 
(58) 



Unfortunately we cannot solve these equations for the general situation (except for the 
trivial case L = 2 and L = 3 - see below). It is much easier to consider non-minimal 
critical points of ||^|| 2 . In order to find states characterizing other SLOCC classes it is 

(59) 



sufficient co consider solution to the problem 

a*v = Xv 



where a* and a — 



A 

-A 



with Ae (0,1]. One easily computes that 



a*\k, L) = Xk\k, L) , fc = 0, 1,...,L, 



(60) 



where A^ = (2k — L) A. For the considered values of A all eigenvalues A^ are non- 
degenerate and therefore the only candidates for critical states are Vk- One checks that 



*([«*]) = 



2k-L 
2L 









L-2k 
2L 



(61) 



Clearly, for k = |~|f] . . . L, the mapping <P ([ufe]) gives a of the desired property. There- 
fore we have the following states parameterizing inequivalent SLOCC classes 



Vk 

Direct computations show that 



\k, L),k = 



(62) 



{T [vk] G.[v k ]Y 



Span{|m, L), i\m, L) , m € {0, 1, . . . , k - 2} U {k + 2, . . . , L}} 

(63) 

where it is assumed that k > 2 and k < L — 2. States v and v\ do not appear as critical 
states with e «t + . One has to consider separately the cases k = L — 1 and 

k = L, i.e. 



(7k]G.M)- 



Span {|to, L), i|m, L) , m e {0, 1, . . . , fc - 2}} , (64) 



for fc = L — 1 and k = L. Remarbly not all directions from (TL^C?.^]) " cause the 
decrease of ||^|| 2 . Closer look at the Hessian of function (M) = W ^ 

at [w] — [vk] restricted to (Tj^jG.fvfc])^" reveals that ||^|| 2 decreases in the directions 
given by vectors \k, L), i\k, L) for k € {0, . . . , [^J }. Therefore, 

~ L~ 



md([v k ]) 



(65) 
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6.7. Three five-state fermions. As the last example let us consider a system consisting 
of three fermions with the five-dimensional single particle Hilbert spaces. It follows 
that U = A 3 (C 5 ), K = SU (5), G = SL(5, C). In C 5 we fix the orthonormal basis 

{|1), |2), |3), |4), |5)} and chose the following (also orthonormal) basis of A 3 (C 5 ) 
(see the notation in Equation ([38))): 

|*i> *3), 1 < k < h < k < 5 (66) 
Any state v has the decomposition 

v = oi|l,2,3) +0211,2,4) +oa|l,2,5) + a 4 |l,3,4) +a 5 |l,3,5)+ (67) 
+a 6 |l, 4, 5) + 6i]2, 3, 4) + 6 2 |2, 3, 5) + b 3 \2, 4, 5) + ci|3, 4, 5) , (68) 

where the scalar coefficients fulfill the normalization condition. In order to find p< v ] one 
has to compute terms of the form trj (|1, 2, 3)(2, 3, 4|) (see Eq. (pTji). When two indices 
in states that appear in such an expression coincide the result is nonzero, 

tri ( |ii , i 2 , * 3 >0'i , J2 , J's I ) ^0 (69) 
t (70) 
the set {ii, i 2 , 13} (1 j'2, J3} has more than one element. (71) 

It is then easy to verify that for each basis vector z 2 , 13) there are six other basis 
vectors that give non vanishing cross terms. For an exemplary |1. 2, 3) these are: 

|1,2,4), |1,2,5), |1,3,4), 11,3,5), |2, 3, 4), |2, 3, 5). (72) 

A system in H. = /\ 3 (C 5 ) can be considered as consisting of two holes and the struc- 
ture of the spectrum of the rescaled density matrix p = 3 • p\ v ] is known |28 1. 

- All eigenvalues of p are bounded above by 1. 

- The largest eigenvalue X max — 1. 

- All other eigenvalues have at least double degeneracy. 

Using these properties it can be easily shown that there are only two nonequivalent 
critical sets of ||/x|| 2 . They are parameterized by the states 

wi = 11,2,3), (73) 

u 2 = -^ (| 1,2,3) + 1 1,4, 5}) . (74) 

In order to compute corresponding indices we check that 

(T K] G.[^]) Xa, =Span{|l,4,5), i|l,4,5), |2,4,5),i|2,4,5), |3,4,5), i|3,4,5)} , 

(75) 

and 

(T [v2] G.[v 2 ]) ±u = 0. (76) 

We conclude that the the Morse index of [v±] is equal 6. On the other hand the orbit 
G.[i> 2 ] has the maximal dimension (8) and thus index of [v?\ vanishes. 
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7. Summary 

A geometric description of entanglement, or more generally, quantum correlations pre- 
sented in our paper proved to be very fruitful in analysis of various aspects of composite 
systems. The main achievement is a universal algorithm allowing for determination of 
inequivalent SLOCC classes of pure states. For qubit systems the algorithm is effective 
- we showed how it works for two-, three-, and four-qubit systems. Since the canonical 
form of a state is known for arbitrary number of qubits, similar calculations can be, in 
principle, performed for such systems, although definitely with more effort. In any case 
the algorithm can be implemented numerically. 

The algorithm distinguishes in a natural way two types of SLOCC classes: those 
which are parameterized by nontrivial closed G-orbits where G is the group of SLOCC 
transformations, and those which belong to the null cone N, i.e. contain in their closures 
as a unique closed G-orbit. The orbits of the first type exhaust almost all G-orbits. 
From this point of view orbits in the null cone seem to be not particularly interesting, 
but in fact it is not so. Indeed, starting from the equivalence classes of pure states one 
can construct a classification of mixed quantum states. Different classes of mixed states 
are then defined as convex sums of projectors on pure states from a particular class. Pure 
states classes of measure zero may give rise an the level of mixed to sets of non-zero 
measure l30ll . 

The main tool used in our approach was the momentum map and its geometric prop- 
erties - in particular convexity of its image (or its appropriate restriction). The momen- 
tum map is a natural construction whenever a group acts on a symplectic manifold in a 
symplectic manner. Such a situation is characteristic not only for the case of quantum 
entanglement for distinguishable particles where the appropriate group consist of all 
local unitary transformations acting the space of multiparticle states, but mutatis mu- 
tandis also in the bosonic and fermionic cases, when additional symmetry conditions 
are imposed upon states. The geometric and group-theoretic origins of our algorithm 
make it versatile enough to cover also these cases. 

In the course of argumentation we exhibited also some additional interesting proper- 
ties of the momentum map useful in characterizing entanglement of states by appropri- 
ately defined distance of the SLOCC orbits from the origin and by counting directions in 
which entanglement can be increased under non-SLOCC operations. Finally we showed 
its connection with potentially measurable physical observables. 
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